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Abstract 

A simple model of nonfemtoscopic particle correlations in proton-proton collisions is proposed. 
The model takes into account correlations induced by the conservation laws as well as correla- 
tions induced by minijets. It reproduces well the two-pion nonfemtoscopic correlations of like-sign 
and unlike-sign pions in proton-proton collision events at yfs = 900 GeV analyzed by the ALICE 
Collaboration. We also argue that similar nonfemtoscopic correlations can appear in the hydrody- 
namic picture with event-by-event fluctuating nonsymmetric initial conditions that are typically 
associated with nonzero higher-order flow harmonics. 



PACS numbers: 13.85.Hd, 25.75.Gz 



I. INTRODUCTION 



The method of identical particle correlation femtoscopy gives us the possibility to measure 
length and time scales with an accuracy of more than 10 -14 m and 10~ 22 sec, respectively. 
It is widely used now for the determination of sizes and lifetimes of the sources of particle 
emission such as the systems created in heavy ion, hadron, and lepton collisions (for reviews 
see, e.g., Ref. [lj]). The method is grounded on the Bose-Einstein or Fermi-Dirac symmetry 
properties of quantum states and, in fact, measures correlations between numbers of identical 
particles with close energy and momentum. 

In the method of correlation femtoscopy the space-time structure of the systems is usually 
represented in terms of the interferometry radii that are the result of a Gaussian fit of a two- 
particle correlation function depending on the momentum difference in the pair. For bosons 
such correlations, reflecting spatiotemporal scales of an extended source, are caused by the 
quantum Bose-Einstein statistics. These correlation functions can be obtained from the 
ratio of the two- (identical) particle momentum spectra to the product of the single particle 
ones. The former typically need to be corrected for the Coulomb and strong final state 
interactions (FSI) that depend on space-time points of the last collisions for the detected pair 
(see, e.g., the recent review [2] and references therein). In pioneering papers jsj] the measured 
interferometry radii were interpreted as geometrical sizes of the systems. Later on it was 
found 4|, |5] that because the typical systems formed in experiments with heavy ion collisions 
are expanding and have inhomogeneous structures, the above geometrical interpretation 
is not complete, and the interpretation of interferometry radii as homogeneity lengths in 
the systems was proposed [6(. It was demonstrated that the homogeneity lengths carry 



information about the dynamical properties 
and depend on the mean momentum of pairs 



rate of expansion, lifetime, etc.) of the source 



IE 



6|. In particular, it is firmly established that 



specific transverse momentum dependence of femtoscopy scales - interferometry radii - in 
heavy ion collisions is mostly caused by a collective (hydrodynamical) transverse expansion 
of the systems formed in these collisions [lj]. 

The situation is more complicated, and the above method has to be modified for elemen- 
tary particle collisions, like p + p, which have smaller spatiotemporal scales as compared to 
heavy ion collisions. It became clear jl-js] that for relatively small systems the additional 
two-particle correlations affect the correlation functions in the kinematic region where quan- 



turn statistical (QS) and FSI correlations are usually observed. The well-known example of 
such additional correlations is the correlation induced by total energy and momentum con- 



servation laws (see, e.g., Ref. [10|)- As opposed to the QS and FSI correlations, which are 
familiar from the correlation femtoscopy method, and so are sometimes called femtoscopy 
correlations, these correlations are not directly related to the spatiotemporal scales of the 
emitter and are therefore called nonfemtoscopic correlations. Since the latter noticeably 
affect correlation functions for small systems, the interferometry radii extracted from the 
complete correlation function in p + p collisions depend strongly on the assumption about 
the so-called correlation baseline - the strength and momentum dependence of the nonfemto- 
scopic correlations jz-9]. It has an influence on the interpretation of the momentum depen- 
dence of the interferometry radii in p + p collisions, where the possibility of hydrodynamic 
behavior of matter is questionable. Therefore, for successful and unambiguous applications 
of the correlation femtoscopy method to elementary particle collisions, one needs to know the 
mechanisms of nonfemtoscopic correlations to separate the femtoscopic and nonfemtoscopic 
correlations. 

Recently, the ALICE Collaboration utilized some event generators, which do not include 
effects of quantum statistics, for an estimate of the correlation baseline (i.e., nonfemtoscopic 







correlation function of identical pions) under the Bose-Einstein peak at LHC energies 
It was motivated by a reasonable agreement of the corresponding event generator simulations 
with the experimental data for correlation functions of oppositely charged pions in p + p 



collisions at the same energy 



8] . The calculated correlation baseline has been utilized by 



the ALICE Collaboration to extract femtoscopic correlations from measured identical pion 



two-particle correlation functions {j]. 

Because the utilized event generators account for energy-momentum conservation and 
emission of minijets, it was conjectured in Refs. {?], [§]] that some specific peculiarities of 
the unlike-sign pion correlations as well as like-sign nonfemtoscopic pion correlations can 
be caused by the jetlike and energy-momentum conservation induced correlations. In what 
follows, we support this conjecture. We develop a simple analytical model with the minimal 
number of parameters for the two-pion correlations induced by minijets and transverse mo- 
mentum conservation law, and show that this model can fit the correlations of unlike-sign 
p,o n p alre a. V. = 900 GeV P + P comsions _ b y tbe ALICE Coloration fl. Also, 
with a reasonable change of parameters, the model can fit nonfemtoscopic correlations of 



like-sign pion pairs obtained from the event generator simulations of p + p collision events 
at the same energy that are reported in Ref. Our model is simple and analytical, and 
clearly demonstrates the interplay between minijet and conservation law induced correlations 
in the formation of the nonfemtoscopic correlations. 

But this is not the whole story. Despite the fact that the event generators do reproduce 
the unlike pion correlations as has been verified by the ALICE Collaboration, there is enough 
room to doubt whether the nonfemtoscopic correlations of like-sign pions are properly sim- 
ulated by them. Indeed, none of the utilized event generators can reproduce the LHC data 



on the multiplicities and momentum spectra well (see, e.g., 12| and references therein). 
This suggests that some essential ingredients may be missing in these event generators. The 
possible candidate for missed dynamics is hydrodynamics. The latter describes well the dy- 
namics of heavy ion collisions, see, e.g., Ref. 



13] and references therein. If hydrodynamics 
is applied also for p + p collisions with high multiplicities [l^] , then the other mechanisms 
of nonfemtoscopic correlations should be taken into account. In what follows, we present 
some heuristic arguments, based on an illustrative analytical model, that unlike-sign two- 
pion correlation functions calculated in hydrodynamics with event-by-event asymmetrically 
fluctuating initial densities can be qualitatively similar at relatively low q inv to the ones 
calculated in PHOJET-like event generators, where these correlations for relatively low qi nv 
are mainly caused by the minijets. Finally, we briefly discuss what these results can mean 
for modeling of the correlation baseline (i.e., identical pion nonfemtoscopic correlations) and 
for applications of the correlation femtoscopy method to p + p collisions. 

II. DEFINITIONS AND PARAMETRIZATIONS OF TWO-PARTICLE CORRE- 
LATIONS 

The two-particle correlation function is defined as 

where P{pi,P2) is the probability of observing two particles with momenta p x and p 2 , while 
P(pi) and P(j>2) designate the single-particle probabilities. Experimentally, the two-particle 



11] simulations 



1 For convenience, we compare results of our model with the PHOJET event generator 
reported in Ref. [7[. Note that the simulations carried out by the ALICE Collaboration gave similar 
results for all utilized event generators Q, 0] . 



correlation function is defined as the ratio of the distribution of particle pairs from the 
same collision event to the distribution of pairs with particles taken from different events. 
In heavy ion collisions almost all the correlations between identical pions with low relative 
momentum are due to quantum statistics and final state interactions. In this case 

CW 2 ) = CV(p,q), (2) 

where p = (pi + P2)/2, q = P2 — Pi, and Cf denotes the femtoscopic correlation function. 
For identical bosons Cf is often parameterized (after corrections for FSI correlations) by 
the Gaussian form which for the one-dimensional parametrization looks like 

C F (\p\,q inv ) = 1 + Aexp (-R 2 iTlv q 2 inv)- (3) 
Here A describes the correlation strength, Ri nv is the Gaussian "invariant" interferometry 



radius, and q inv = (p 2 — pi) 2 — (E 2 — Ei) 2 is equal to the modulus of the 3-momentum 
difference in the pair rest frame. 

In elementary particle collisions additional (nonfemtoscopic) correlations, like those aris- 
ing from let /string fragmentation and from energy and momentum conservation (see, e.g., 



Refs. 

erty, 



can also give a significant contribution. Then, assuming the factorization prop- 



C(pi,p 2 ) = CV(p,q)CW(p,q). (4) 

Here Cnf denotes the nonfemtoscopic correlation function, and in the simplest case the 
nonfemtoscopic effects can be parametrized as, e.g., 2nd order polynomial, 

Cnf ( I P I , qinv) =a + bq inv + cqj mv . (5) 

This form can be used together with some parametrization of Cp (e.g., with dSJ)) in order to 
fit the correlation function C(pi,p2) for small systems, as has been done, for example, by the 
STAR Collaboration for two-pion correlation functions in p + p collisions at \fs = 200 GeV 
^ . At c > the phenomenological parametrization (jHJ) explicitly reproduces the well known 
effect of positive correlations between particles with large relative momenta |q| caused by the 
energy-momentum conservation laws; see the energy and momentum conservation- induced 
correlation model for C^f |lQj- Note that a, b, and c in Eq. flS} depend, in general, on |p|. 



III. PROBABILITY DENSITIES OF DISTINGUISHABLE EQUIVALENT PAR- 
TICLES 



Aiming to estimate the nonfemtoscopic pion correlations, we consider emitted pions as 
distinguishable, yet equivalent noninteracting particles with symmetrical probability density 
functions, thereby, excluding femtoscopic QS and FSI correlations. To a certain extent this 
corresponds to the (quasi) classical approximation used in current event generators, like 
PHOJET. 

Let us assume that N bosons of the same species (say, pions) are produced with momenta 
p x , p N in (N + X) multiparticle production events. Then the iV-particle probability den- 
sity P/vG°i ; ■■■iPn) is a symmetrical function for all N\ permutations of the particle momenta 
Pi. For convenience, it is normalized by 

J dQ p P N (px, ...,p N ) = 1, (6) 
where dVt p = ^^...^ IL . Then the iV'-pion probability density, N' < N, is defined as 

P N ( Pl , ...,p N ,) = j <m p *E*5®{ Pl - p* 1 )...E* N ,6®(p N , - P* N ,)MP*n -,P*n)- (7) 

In what follows, we use an assumption of distinguishability of equivalent particles which 
means the absence of quantum interference between possibilities corresponding to all AH 
permutations of the particle momenta p^. Therefore the symmetrized iV-particle probability 
density is defined as 

l - r 

P N ( Pl , ..., PN ) = — / dn p *E*5 {3) ( Pl - p*)...E* k 5^(p N - p* k )P N (pl ...,p* N ), (8) 

where the nonsymmetrized iV-particle probability density Pn(pi, ■■-,Pn) is normalized to 
unity, and A^! in the denominator is required to guarantee the normalization condition ([6]). 
Taking into account Eq. (J7|), we see that the single-particle probability Pjy{pi) and the 
two-particle probability Pn{vuP2) can be written as 

Pn(Pi) = E / «V^^(Pi - Pistil ...,P* n ), (9) 
i=l J 

i - r 

Pn(Pi,P2) = N{N _ 1} E J dn p *E*6®( Pl - p*)E;5^(p 2 - p*)P N (pl ...,p* N ). (10) 



The nonsymmetrized iV-pion probability density in such events reads 

If NX 

P n {pi,-,Pn) = x22 / <MkS i %a+Pb-^Pi-22k j )\M N+x (j> 1 ,...,k x )\\ (11) 

X J i=l j=l 

where M N+X (pi, k x ) is the nonsymmetrized (N + X)-particle production amplitude, p a 
and pi, are 4- momenta of colliding particles (protons in p + p collision events), and K is the 
normalization factor, 

» NX 

k = J2 dn k dn p 6^\p a + p b -J2^-J2 k ^ MN + x ^---' kx ^ 2 - ( 12 ) 

X J i=l j=l 

Expression ( ITT]) for Pn(Pi, ■■■,Pn) is rather complicated because, in particular, it depends 
on X particles that are produced in addition to N pions. Production of additional parti- 
cles also means that one can hardly expect that the total energy or momentum of the pion 
subsystem are constants in the system's center of mass; instead, one can expect that they 
fluctuate from event to event. Then, motivated by Eq. (jlip . we assume that a nonsym- 
metrized iV-pion probability density can be written as 

P N (px, ...,p N ) = ^(Pi» -;Pn)Fn(Pi, -,Pn), (13) 

where F N (p 1: ...,p x ) is a nonsymmetrized function of pionic momenta and 5(pi, ...,Pn) de- 
notes the average constraints on the iV-pion states that appear due to energy and momentum 
conservation in multiparticle production events. Then the normalization factor is 

K = J d£L p 6(pi,...,p N )F N (pi,:.,PN)- (14) 

If the only correlations are the correlations associated with energy and momentum conser- 
vation, we have 

F N (p u ...,p N ) = f(px)f(p2)-f(PN-i)f(pN), (15) 

and calculations in the large N limit of single-particle and two-particle probability densities 
result in the special case of the EMCIC parametrization [10] of the correlations induced by 
the energy and momentum conservation laws. However, such a simple prescription cannot 
result in the unlike-sign pion correlations measured by the ALICE Collaboration 7[ • Also it 
does not reproduce the nonfemtoscopic like-sign pion correlations generated by the PHO JET 
event generator 



IV. ANALYTICAL MODEL FOR THE TWO-PION CORRELATIONS INDUCED 
BY MINUETS AND MOMENTUM CONSERVATION 



[7 



The ALICE Collaboration has analyzed the unlike-sign pion correlations |7| and found 
that they can be well reproduced by event generators which account for, among others 
factors, total energy-momentum conservation and minijet production. In what follows, we 
assume that just these two factors induce the observed behavior of the unlike-sign pion 
correlations and are responsible for the nonfemtoscopic correlations in like-sign ones. In our 
simple model we start from zero total transverse momentum of the system, keeping in mind 
that for a subsystem this statement should be weakened. Also, we neglect the constraints 
conditioned by the conservation of energy and longitudinal momentum, supposing that the 
system under consideration is an iV-pion subsystem in a small midrapidity region of the 
total system. Then in the first approximation 

6{pi, ...,p N ) = 5 (2) (Pn + p T2 + ... + Ptn), (16) 

where Pti, Pt2, ■■■Ptn are transverse components of the momenta of the N particles. 

Let us assume that there are no other correlations in the production of iV-pion states 
except the correlations induced by the transverse momentum conservation and cluster (mini- 
jet) structures in momentum space. For the sake of simplicity we assume here that only the 
two-particle clusters appear. Then one can write for fairly large N ^> 1 

F N (pi, ...,p N ) = f(pi)-f(pN)Q(pi,P2)-Q(pN-i,PN), (17) 

where Q(j>i,Pj) denotes the jetlike correlations between momenta p { and p,,; the existence 
of such correlations means that Fn cannot be expressed as a product of one-particle dis- 
tributions. Then, utilizing the integral representation of the ^-function by means of the 
Fourier transformation, 5^(pt) = (2rr)~ 2 J d 2 r T exp(ir T p T ), and accounting for Eqs. (j5J), 
(|T3|) . (fTB|) . (JTTj) . the single-particle probability reads 

Pn{pi) = 7>W / d2rT ° N ^ rT ^ ( 18 ) 

where 

G N (p l ,r T ) = J d^^^ 3 )(p 1 -pDe irT ^ + --- +p ^F JV (K,...,^). (19) 



A possibility of different cluster configurations of particles means, in particular, that reg- 
istered particles with momenta p x and p 2 can belong either to different minijets or to the 
same minijet. Then, taking into account Eqs. (TTTJj) . ( 1T5]) . (TIT)]) . ( TT71) . we get 

N uV«u. „ , , N(N - 1) - N p2jet/ 



where 



Pn(PuP2) = n{n _ 1) P n(PuP2) + V(iV-l) P » ^ 

PN Jet (PuP2) = I d%G^( Pl ,p 2) r T ), (21) 

P 2 J et (puP2) = J rf 2 r T Gr(Pi,P2,r T ), (22) 



and 



(23) 



G# et ( Pl ,p 2 ,r T ) = J dn p *E*5W( Pl - pl)E*5^(p 2 - V * 2 )e^^-^T N ) p Ni 

Gf i ( P i,p 2 ,r T ) = J ^E* ( 5( 3 Hp 1 -pDE; ( 5( 3 )(p 2 -p;)e^^ + --- +p ^F w . (24) 

Here F N = F N (pl, ...,p* N ). The first term on the right-hand side of Eq. f j20l) is associated 
with events where the two registered particles belong to the same minijet, and the second 
term corresponds to events where the particles are from different minijets. Evidently, the 
former is relatively rare; however, notice that the first term can be significant for small 
systems with not very large N. 

Now let us check whether this model can reproduce, with reasonable parameters, the 
correlation functions of unlike-sign pions measured by the ALICE Collaboration [3] and 
nonfemtoscopic correlations of like-sign pions that are generated in the PHO JET simulations 
and utilized as the correlation baseline by the ALICE Collaboration [7]. Calculations within 
the model will deliberately be as simple as possible just to demonstrate its viability Here 
we do not use the approximate methods like the saddle point approach; instead, we utilize 
appropriate analytical parametrizations of the functions of interest, namely, 

f( Pi ) = Et exp (-^r) ex P {~Jf) ' ( 25 ) 

and 

Q( Pi ,p j )=exp(- {Pi -J j)2 \ (26) 

where T T , T L , and a are some parameters, and in what follows, we assume that T L ^> T T . 
In accordance with the ALICE baseline obtained from the PHOJET event generator simu- 
lations, we assume that only q inv is measured for each p T bin. Assuming that longitudinal 



components of the registered particles are equal to zero, pn = P2L = 0, we approximate qf nv 
as 

, , / m 2 + pg,sin 2 (/. \ 

where denotes the unregistered angle between p T and q^, Prqr = IptIIQtI cos0. Then 

d(f)P N (p 1 ,p2 
J 27r d(f)P N (p 1 )P N (p 2 

and, taking into account Eq. (1201) . we get 

(29) 



^JVi^ilPrl) Qinv) — r2 7r ; , ~ / — N ~ ; — 7; l/°J 



where 



Cr( |p r |, fc . ) = Xi^%^., (30) 



2 



CF(M.*-)= {f d0Pr(p " ft) . (31) 

/„ d(t>P N {pi)P N {p 2 ) 

It is well known (see, e.g., Ref. 10]) that the influence of exact conservation laws on 
single-particle and two-particle momentum probability densities at the iV-particle produc- 
tion process depends on the value of N and disappears at N — > oo. Since one considers 
a subsystem of iV pions but not the total system, to weaken the influence of the total 
transverse momentum conservation on pions we shall consider C]^ et and C^ et with M > N 
instead of C]^ et and C^ et in Eq. (1291) . This is the simplest way to account for a weakened 
conservation law in our model. At the same time, the factor 1/(N — 2) in ( )29|) remains the 
same since it is associated with the combinatorics of the distribution of particles between 
clusters in momentum space ("minijets"), which happens whether or not one weakens the 
total momentum conservation law. Also, for more exact fitting of the data points in each 
average transverse momentum bin, we utilize the auxiliary factors A. Then Eq. ( 1291) gets 
the form 

= A(|p T |)(Ctf(|p T |,g_) + — C^Wvrlq^)). (32) 

The results of our calculations of the nonfemtoscopic correlation functions CV p are shown 
in Figs. [THTUlin comparison with correlation functions reported by the ALICE Collaboration 
[?| for different transverse momenta of pion pairs (actually, we performed calculations for 

10 



the mean value in each bin). The auxiliary factors A(|pr|) differ from unity only slightly^ 
The data for unlike-sign pion correlations measured by the ALICE Collaboration as well as 
for the PHOJET simulations of like-sign two-pion nonfemtoscopic correlation functions at 
midrapidity for the total charged multiplicity N c h > 12 bin in p + p collisions at y/s = 900 
GeV are taken from Refs. {?| and 15|. Note that correlations of nonidentical pions measured 



by the ALICE Collaboration, as well as the PHOJET simulations of identical two-pion 
correlation functions, demonstrate Coulomb FSI correlations at the lowest qi nv bin and 
peaks coming from resonance decays. These Coulomb FSI, as well as contributions from 
resonance production, are not taken into account and so are not reproduced in our model. 

The presented results are obtained for M = 50, Tt = a = 0.65 GeV (to minimize the 
number of fit parameters, we fixed Tt = a for all calculations)!^, and the fitted values of A" 
are different for like-sign and unlike-sign pion pairs, namely, N ±ziz = 20 for the former and 
A rH =11 for the latter. The relatively high value of M can be interpreted as a residual effect 
on the pion subsystem of total energy-momentum conservation in a multiparticle production 
process. The relation A^ + ~ < N ±i: between fitted A" values means that the magnitude of the 
correlations induced by a minijet for unlike-sign pion pairs is higher than for like-sign ones. 
This happens because in the former there is no local charge conservation constraint for the 
production of oppositely charged pion pairs and, therefore, one can expect less identically 
charged pion pairs from the fragmenting minijets than oppositely charged ones. 

One can see from the figures that the behavior of the nonfemtoscopic correlation functions 
of pions, Cnf, is reproduced well despite the simplicity of our model. This is a result of the 
competition of the two trends: an increase of the correlation function with qi nv because of 
momentum conservation and a decrease of it due to fragmentation of one minijet into the 
registered pion pair. Figures [5] and (TU] also demonstrate the relative contribution of the first 
and second terms in Eq. (1321) to the nonfemtoscopic correlation functions. 



2 Namely, they are 0.95, 0.96, 0.98, 0.99, 0.99 for like-sign pions (in Figs. HE respectively), and 0.92, 0.92, 

0.93, 0.91, 0.82 for unlike-sign pions (in Figs. IMTUI respectively). 

3 Note that with these parameter values the mean transverse momentum (px) is about 0.58 GeV. 
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V. ANALYTICAL MODEL FOR EVENT-BY-EVENT MOMENTUM SPECTRA 
FLUCTUATIONS 



As it follows from our previous discussion, the lower magnitude of the nonfemtoscopic 
correlations at relatively low qi nv for like-sign pion pairs as compared to the correlations of 
unlike-sign pions is natural if pions are produced through minijet fragmentation. However, 
this cannot be the case for other production mechanisms that do not include noticeable 
production of minijets. For example, if thermalization takes place in p + p collisions and 
hydrodynamical evolution forms the "soft" momentum spectrau, then the production of mini- 
jets at relatively low px is, typically, reduced. In this case utilization of the nonfemtoscopic 
correlations of like-sign pion pairs obtained in the PHO JET and similar event generators as 
a correlation baseline for the femtoscopic correlations, see Eq. (HI) , can be in doubt. Then 
the question arises whether the unlike two-pion correlations in p + p collisions, which are re- 
produced in PHOJET-like models, are ultimately caused by minijets and conservation laws 
only, or whether a similar behavior can be attributed to hydrodynamics also. 

First, note that there are no correlations induced by the exact global energy-momentum 
conservation in hydrodynamic models, and corresponding conservation laws are satisfied 
only in average for particles that are produced at some hypersurface where hydrodynamics 
is switched off. The global energy-momentum conservation constraints can be added on an 
event-by-event basis if the post-hydrodynamical hadronic stage is calculated by means of 
some hadronic cascade model (the so-called "hybrid" model). One more source of the non- 
femtoscopic correlations in such models are event-by-event fluctuations of initial conditions 
for the hydrodynamical stage. These fluctuations result in fluctuations of the two-particle 
and single-particle momentum spectra, and, as usual, the effect of the fluctuations is more 
pronounced for small systems. Then, an important question is whether such correlations 
can be similar to minijet induced correlations that, as we know, may reproduce unlike-sign 
pion correlations at relatively low qi nv . 

Let us give an illustrative example of nonfemtoscopic correlations that appear because of 
event-by-event fluctuating initial conditions for the hydrodynamic stage (in hybrid models 
this stage is matched with the subsequent hadronic cascade stage) and are similar to the 



4 Hydrodynamic models, perhaps, can give a reasonable description of elementary particle collisions; see, 

e.g., y . 
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ones produced by minijets. Suppose that the iV-particle probability density is defined as 

Pn(puP2, -,Pn) = / J w(u i )P N {pi,p 2 , ...,p N ;Ui), (33) 

i 

where Pjv(pi,P2, ■■■,Pn', u%) is the iV-particle probability density for some u% type of the 
initial conditions, and w(ui) denotes the distribution over initial conditions, Y^i w ( u i) = 1- 
To analyze the possible effect of fluctuating initial conditions, here we neglect conservation 
law constraints and the production of minijets. Because we assume uncorrelated particle 
emissions for each specific initial condition, one can write 

P N (Pi,P2, ...,p N ;Ui) = f(p 1 ;u i )f(p 2 ;Ui)...f(p N ^ 1 ;u i )f(p N ;u i ), (34) 

where we normalize f(p;Ui) as J ^/(p;«i) = 1, and then K — 1; see Eqs. ([TBI . (1T4"]) . The 
two-particle nonfemtoscopic correlation function Cnf then reads 

C NF {p u p 2 ) = ^ — r — TT r — 77 r. (35) 

Li ™{ui)f{Pi; Ui) w(uj)f(p 2 ; Uj) 

Evidently, the different type of fluctuation, i.e., the form of the distribution w(ui), leads 
to a different behavior of the nonfemtoscopic correlations. To illustrate that fluctuations 
can lead to the nonfemtoscopic correlation functions that are similar to the ones induced by 
minijets, let us consider the toy model where 

2 

w(u T ) = — exp(— u^a 2 ), (36) 

71 

f(p;u T ) = ^Eexp(-(p T - UT ) 2 /3 2 )exp(-p| 7 2 ), (37) 



7T 



and normalization is chosen in such a way that J c1 2 utw(ut) = 1 and J ^/(p; i*r) = 1. 
The main feature of such a model is that event-by-event single-particle transverse momen- 
tum spectra have a maximum for event-by-event fluctuating p T values. Such momentum 
spectrum fluctuations could take place, e.g., in hydrodynamics with a highly inhomogeneous 
initial energy density profile without cylindrical or elliptic symmetryjj One can easily see 
that in this case CVf decreases with gf,, 

/3 4 

C NF (p, q) ~ ex P(~ 2 ^ 2 + ; g2) gr), (38) 



5 It seems that this is the case in heavy ion collisions, where nonsymmetrical fluctuations of initial conditions 
lead to nonzero v$ and higher flow harmonics (see, e.g., Ref. |17j). 
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and this means (after taking into account ff271) and ([28]) ) that Cnf decreases with qf nv too, 
which is similar to the behavior of Cnf if the nonfemtoscopic correlations are induced by 
minijets. At the same time, unlike the latter, the hydrodynamical fluctuations lead to similar 
correlations for like-sign and unlike-sign pion pairs. Then, our analysis suggests that, up to 
different resonance yields, the value of the slope of the correlation baseline at relatively low 
qi nv can be somewhere between pure hydrodynamic (i.e., the same as for nonidentical pion 
pairs) and pure minijet (i.e., lower than for nonidentical pion pairs) scenarios. 



VI. CONCLUSIONS 

We presented here a simple analytical model that takes into account correlations in- 
duced by the total transverse momentum conservation as well as correlations induced by 
the minijets. It is shown that the model gives a reasonable description of the correlations of 
nonidentical pions measured by the ALICE Collaboration [7] in p + p collisions at yfs = 900 
GeV, and also the nonfemtoscopic correlations of identical pions generated in the PHOJET 
simulations of p + p collisions and utilized by the ALICE Collaboration [3] as the correlation 
baseline. We conclude that the cluster (minijet) structures in the final momentum space 
of produced particles can result in noticeable nonfemtoscopic two-pion correlation functions 
that decrease when qi nv grows at relatively low qi nv , while the global energy- momentum con- 
servation constraints typically result in an increase with qi nv for fairly high qi nv . Our model 
can be utilized for simple estimates of the nonfemtoscopic correlations induced by minijets 
and conservation laws that contribute to the total two-particle correlation functions. 

There can be different types of multiparticle production mechanisms, and some of them 
could result in qualitatively similar nonfemtoscopic correlation functions. We presented 
heuristic arguments that the two-pion nonfemtoscopic correlation functions calculated in 
hydrodynamics with event-by-event fluctuating initial conditions can be qualitatively simi- 
lar at relatively low q inv to the ones calculated in the PHOJET-like generators, where the 
nonfemtoscopic correlations for low qi nv are mainly caused by minijets. It is worth noting 
an important difference between the nonfemtoscopic correlations induced by minijets and 
hydrodynamical fluctuations: while the former lead to a higher magnitude of the nonfemto- 
scopic correlations for unlike-sign pion pairs as compared to like-sign pions, the latter result 
in a similar (up to the resonance contributions) strength of the nonfemtoscopic correlations 
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for identical and nonidentical pions. Then, if the applicability of hydrodynamics to p + p 
collisions is supported, such an analysis allows one to estimate the correlation baseline and, 
so, to extract the femtoscopic scales in these collisions by means of tuning the hydrokinetic 
model to reproduce the experimental unlike-sign pion correlations. 

Because the particle production mechanisms in p + p are still unclear, the nature of 
nonfemtoscopic correlations in these collisions is also an open question. Different dynamical 
models, that reproduce unlike pion correlations, can give different estimates of the correlation 
baseline and, so, lead to different results for the correlation femtoscopy analysis of the space- 
time scales of the collision process. In our opinion, this difficulty can be overcome if the 
correlation analysis is applied not to femtoscopic correlation functions, but to experimental 
data reported for the total correlation functions of pion pairs. Then the space-time scales 
can be estimated by means of the dynamic model that will also be able to reproduce, among 
other observables, these compete correlations. 
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FIG. 1. The nonfemtoscopic correlation functions of like-sign pions in the 0.1 < px < 0.25 GeV 
bin from a simulation using PHO JET B Q (solid dots) and those calculated from the analytical 



model: minijets + momentum conservation (solid line). See the text for details. 




FIG. 2. The same as Fig. []]but in the 0.25 < p T < 0.4 GeV bin. 
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FIG. 3. The same as Fig. [Qbut in the 0.4 < p T < 0.55 GeV bin. 
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FIG. 5. The nonfemtoscopic correlation functions of like-sign pions in the 0.7 < px < 10 GeV 
bin from a simulation using PHOJET (solid dots) and those calculated from the analytical 

model (solid line) . The contributions to the nonfemtoscopic correlation function from the first term 
of Eq. (j32|) (dotted line) and from the second one (dashed line) are also presented. 




FIG. 6. The correlation functions of unlike-sign pions in the 0.1 < pr < 0.25 GeV bin measured by 
the ALICE Collaboration from Refs. (solid dots) and those calculated from the analytical 

model: minijets + momentum conservation (solid line). See the text for details. 
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FIG. 7. The same as Fig. ©but in the 0.25 < p T < 0.4 GeV bin. 





FIG. 9. The same as Fig. ©but in the 0.55 < p T < 0.7 GeV bin. 
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FIG. 10. The correlation functions of unlike-sign pions in the 0.7 < pt < 1-0 GeV bin measured by 
the ALICE Collaboration from Refs. (solid dots) and those calculated from the analytical 

model (solid line) . The contributions to the nonfemtoscopic correlation function from the first term 
of Eq. ([52]) (dotted line) and from the second one (dashed line) are also presented. 
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